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Abstract 

We consider Abrikosov-type vortex lattice solutions of the Ginzburg-Landau equations of supercon- 
ductivity, consisting of single vortices, for magnetic fields below but close to the second critical magnetic 
field H C 2 — k 2 and for superconductors filling the entire R 2 . Here k is the Ginzburg-Landau parameter. 
The lattice shape, parameterized by r, is allowed to be arbitrary (not just triangular or rectangular). 
Within the context of the time-dependent Ginzburg-Landau equations, called the Gorkov-Eliashberg- 
Schmidt equations, we prove that such lattices are asymptotically stable under gauge periodic pertur- 
bations for k 2 > |(1 — ^-fj) and unstable for k 2 < |(1 — W^y), where /3(r) is the Abrikosov constant 
depending on the lattice shape r. This result goes against the common belief among physicists and 
mathematicians that Abrikosov-type vortex lattice solutions are stable only for triangular lattices and 
k 2 > §. (There is no real contradiction though as we consider very special perturbations.) 

1 Introduction 

1.1 Background 

Macroscopic theory of superconductivity, by now a classical theory presented in any book on superconduc- 
tivity and solid state or condensed matter physics (see e.g. [3TJ[20]), was developed by Ginzburg and Landau 
along the lines of the Landau theory of the second order phase transitions and before the microscopic theory 
was discoverd. At the foundation of this theory lie the Ginzburg-Landau equations for the order parameter 
and magnetic potential. The time-dependent generalization of these equations was proposed by Schmidt 
([28]) and Gorkov and Eliashberg ([13]) and are known as the Gorkov-Eliashberg or Gorkov-Eliashberg- 
Schmidt equations (as well as the time-dependent Ginzburg-Landau equations). The latter equations have 
much narrower range of applicability than the Ginzburg-Landau equations (|31j) and many refinements of 
them were proposed, but even a slight improvement of these equations is, at least notationally, extremely 
cumbersome. 

By far, the most important and celebrated solutions of the Ginzburg-Landau equations of superconduc- 
tivity are vortices and vortex lattices, discovered by Abrikosov ([1]), and known as Abrikosov (vortex) lattice 
solutions. Among other things understanding these solutions is important for maintaining the superconduct- 
ing current in Type II superconductors, i.e. for n > 4=, in our units. 

Abrikosov lattice solutions are extensively studied in physics literature. Among many rigorous results, we 
mention that the existence of these solutions was proven rigorously in [25] [8j [TTJ [6j [32] ((25] [8] [5] deal with 
triangular and rectangular lattice, while [111 151 152" ] . with lattices of arbitrary shape). Moreover, important 
and fairly detailed results on asymptotic behaviour of solutions, for k — > oo and the applied magnetic fields, 
h, satisfying h < ^ logK + const (the London limit), were obtained in [B] (see this paper and the book [37] 
for references to earlier works). Further extensions to the Ginzburg-Landau equations for anisotropic and 
high temperature superconductors can be found in [2] [3] . 

In this paper we address the problem of stability of the Abrikosov vortex lattice solutions within the 
framework of the time-dependent Ginzburg-Landau equations - the Gorkov-Eliashberg-Schmidt equations. 
We consider such solutions for lattices of arbitrary shape (in the extensive literature on the subject such 
solutions are considered only for triangular or rectangular) and for magnetic fields smaller than but close 
to the second critical magnetic field H C 2 — k 2 (the other case of magnetic fields larger than but close to 
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the first critical magnetic field H c \ treated in the literature is not addressed here) and, as common, for 
superconductors filling the entire R 2 . 

We consider the simplest perturbations having the same (gauge-) periodicity as the underlying (station- 
ary) Abrikosov vortex lattice solutions (we call such perturbations gauge-periodic) and prove for a lattice of 
arbitrary shape that, under gauge-periodic perturbations, 

(i) Abrikosov vortex lattice solutions are asymptotically stable for n 2 > |(1 — 

(ii) Abrikosov vortex lattice solutions are unstable for k 2 < |(1 — ~ph))' 

Here /3(r) is the Abrikosov constant depending on the lattice shape t (see Subsection 11.51 for 
the definition). This result belies the common belief among physicists and mathematicians that Abrikosov 
- type vortex lattice solutions are stable only for triangular lattices and n 2 > ^. Here r is complex number 
parametrizing the lattice shapes. (For the definitions of various stability notions see Subsection 1 1.71 ) Gauge - 
periodic perturbations are not common type of perturbations occurring in superconductivity, but the methods 
we develop are fairly robust and can be extended - at the expense of significantly more technicalities - to 
substantially wider class of perturbation, which will be done elsewhere. Moreover, the same techniques could 
be used in other problems of pattern formation, which are ubiquitous in applications. 

To our knowledge, previously, the only known result on stability of the Abrikosov vortex lattice solutions 
concerned orbital stability (under the same type of perturbations) which, though not stated explicitly, can 
be deduced from the variational proof of [25] that the single vortex Abrikosov lattices for k 2 > | are global 
minimizers of the Ginzburg-Landau energy functional on the fundamental cell (see also both authors 
considered only triangular and rectangular lattices; the argument proving the orbital stability for global 
minimizers was proposed in [9]). However variational techniques do not give the asymptotic stability. Also 
variational techniques do not give even the orbital stability for k 2 < | as in this case Abrikosov lattices are 
not global minimizers of the Ginzburg-Landau energy functional on the fundamental cell. 

In the rest of this section we present the basic equations involved, discuss their properties and related 
definitions and present our result. 



1.2 Gorkov-Eliashberg-Schmidt equations 

Macroscopically, states of the superconductors are described by the triples, ('J, A, $) : R + xR 2 — > Cxl 2 xK, 
where "J is the (complex-valued) order parameter, A is the vector potential, and $ is the scalar potential. 
Physically, l^l 2 gives the (local) density of electrons having formed Cooper pairs, curl A is the magnetic field, 
and — dt A — V$ is the electric field. The dynamics is given by the Gorkov-Eliashberg-Schmidt equations (or 
time-dependent Ginzburg-Landau equations) in K 2 which can be written as 

7 a t $* = a j4 * + k 2 (i-|*| 2 )*, 

curl* curl A = -ad t <s>A + Im(§VA*). 

Here n is a positive constant, 7 a complex number with Re 7 > 0, a a positive 2x2 matrix and V4 = V — iA 
and Aa = • V4 are the covariant gradient and Laplacian, and 9t$ is the covariant time derivative 
dt^ = (d t + i<f>)y or d t ^A = d t A + iV<t>, curl A := d Xl A 2 ~d X2 A 1 and curl* / = (d x J,-d x J). The second 
equation is Ampere's law with Jn = — o"(<9 t A+V$) (using Ohm's law with a as the conductivity tensor) being 
the normal current associated to the electrons not having formed Cooper pairs, and J5 = Im(^ , \74^') being 
the supercurrent associated to the electrons having formed such pairs. Note that for a solution (^, A, $) of 
(JT|) the pair (^, A) determines $ through the equation 

A$ = -d t d\vA + divcr _1 [Im(*VA*) - curl* curl A]. (2) 

The equations ([1]) have the structure of a gradient-flow equation for the Ginzburg-Landau energy func- 
tional given by 

£ n (% A) = \j a ||V4*| 2 + I curl A| 2 + y (1 - |*| 2 ) 2 } , (3) 
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where Q is any domain in R 2 . Indeed, (ignoring for the moment the boundary terms) they can be put in the 
form 

d t9 (9,A) = -XE' a ^,A), (4) 

where A is the block-diagonal matrix given by A := diag(7 _1 , c -1 ), and £'q is given by 

£n(*,A) = (-A a *-k 2 (1- |*| 2 )*,curl*curlyl-Im(*V4*)), (5) 

and is formally the L 2 -gradient of £n(^f, A). Though (Q} is not a standard form of the gradient flow, we 
show in Lemma [T2l below that the energy (|3]) decreases under the flow. 



1.3 Ginzburg-Landau equations 

The static solutions of the Gorkov-Eliashberg-Schmidt equations ([1]) are triples (^, A, 0) independent of time. 
In this case (VP, A) satisfies the well-known Ginzburg-Landau equations, which describe superconductors in 
thermodynamic equilibrium and which are given by 

[curl* curlA = Im(Wi*). 
Not surprisingly, they are the Euler-Lagrange equations for the energy functional A) = 0. 

1.4 Symmetries 

The Gorkov-Eliashberg-Schmidt equations ([T]) admit several continuous symmetries, that is, transformations 
which map solutions to solutions: 

Gauge symmetry: for any sufficiently regular function w : R 2 — > R, 

($>(x,t), A(x,t), $(x,i)) ^ (e*"^*(a;,t), A(x,t) + Vn{x,t), $(x,t) - d t r]{x,t)); 
Translation symmetry: for any fcgl 2 , 

{^(x,t), A(x,t), $(£,*)) ^ (^(x + h,t), A(x + h,t), $(x + h,t)); 
Rotation and reflection symmetry: for any R £ 0(2) (including the reflections /(x) 1— > f(—x)) 

(*(x,i), A(x,t), ${x,t)) h-> ($>(Rx,t), R~ l A(Rx,t), <S>(Rx,t)). 

These symmetries restrict to symmetries of the Ginzburg-Landau equations by considering time-independent 
transformations. 

1.5 Abrikosov lattices 

The Abrikosov (vortex) lattice solutions are solutions of the Ginzburg-Landau equations © which exhibit 
double-periodicity in their physical properties. To give a rigorous definition, let £ be a lattice. We say that 
('J, A) is gauge-periodic (with respect to the lattice £), if there exist functions g v : R 2 — > R, v belonging to 
a basis of £, such that for all basis vectors v, 

\y(x + v) =e i9 »Wy{x), _ 



[A(x + u) = A(x) + Vg v {x). 



(In terminology of |29j the pair (ty,A) is equivariant under the group of lattice translations.) An important 
property of gauge-periodic pairs is the quantization of magnetic flux. Let fi be any fundamental cell of C. 
Then the magnetic flux quantization property states that curl A = Inn for some integer n. This can be 
written in terms of the average magnetic flux, b — (curl A)c as b = ^3*, where |fi| denotes the area of fl, 
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and, for / any /^-periodic function, (f)c = Jq /> denotes the average per lattice cell. Using the reflection 
symmetry of the problem, one can easily check that we can always assume n > 0. 

Now any lattice C can be given a basis {i/i,^} such that the complex number r = |^|e* e , where 9 is 
the angle between v\ and ^2, satisfies the conditions that |r| > 1, Imr > 0, — \ < Rer < i, and Rer > 0, 
if |r| = 1. Although the basis is not unique, the value of r is, and we will call it the shape parameter of the 
lattice. 

Note that r, b, and n determine the lattice C up to rotation (which is a symmetry of the Ginzburg-Landau 
equations). We will say that a pair (\I>, A) is of type (r, 6, n), if the underlying lattice has shape parameter r, 
the average magnetic flux per lattice cell is equal to b, and there are n quanta of magnetic flux per lattice cell. 
We also restrict ourselves to C°° pairs (^,A) which suffices for us due to elliptic and parabolic regularity. 
We have the following existence theorem (see |H1 HU 122] )• 

Theorem 1. Let t be any lattice shape parameter and let b be such that k 2 — b is sufficiently small. Then 
there exists an Abrikosov lattice solution = (vl/J, A^) of type (r, b, 1). 

More detailed properties of these solutions are given in Section below. As we deal only with the case 
n=l, we now assume that this is so and drop n from the notation. Note that in this case the average 
magnetic field, b, and the fundamental cell area, |fi|, are related as 

Using the symmetries of the Ginzburg-Landau equations, one can show (see [32] for example) that any 
Abrikosov lattice (^,A) of type (r, b) is gauge equivalent to one satisfying the following conditions: 

(I) (^,A) is gauge-periodic with respect to the lattice CI spanned by r£(l,0) and r£ (Rer, Imr), where 

r _ / 2tt 
'b — Y blmr' 

(II) The function g v (x) in ([7]) can be chosen as g u {x) = \v • Jx, where J is the symplectic matrix 



J 



-1 

1 



and therefore ^ and A satisfy the quasiperiodic boundary conditions 

^{x + v) = e%"- Jx -$(x), and A(x + v) = A{x) + ^ Jv, (9) 

for any element v of a basis of C 
(III) div A(x) = and (A(x) - A° b {x)) C i = °' where A b( x ) ■= l Jx - 

We note that the only continuous symmetry that preserves these properties is the U(l) symmetry 
T 7 : u = (<f, A) T 7 u = (e i7 *, A), V7 G R. 

1.6 Gauge periodic perturbations 

We consider perturbations v = (£, a) of the (r, b)- Abrikosov lattice s.t. + v is again of the type (r, 6). 
Hence it suffices to restrict the problem to any fundamental cell of the lattice • From the definition of 
the type (r, b) pairs, it follows that £ and a satisfy the quasiperiodic boundary conditions 

£(x + v) = e^ l " Jx ^(x) 7 and a(x + u) = a(x), (10) 
for any element v of a basis of C, and with a being divergence- free, and having mean zero, i.e., 



d\v a{x) = 0, and (a(x))ci = 0. 



(11) 
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We introduce our space of perturbations VI to consist of all pairs v — (£, a) 6 i/ 1 (ilJ;C x R 2 ) satisfying 
(fit)]) and (fTTj) . This space is naturally a (real) Hilbert space with the H 1 inner product of v — (£, a) and 
u' = (£', a') given by 

j /• 2 

(«,«V = ToH Re / (&' + V x o^V A o£' + a V + ^Va fe -Va' fe ). 

Here we use that the covariant gradient preserves the quasiperiodic boundary conditions and the dot product 
is in R 2 . We also introduce the L 2 inner product 

(v,v') L 2 = j± r Jte f (tt' + a-a). (12) 

I h I J S~2£ 

Note that (i) if u = ( V I / , A) satisfies (I) - (III), then u — u l satisfies ((TOj) - {EE]); and (ii) the conditions dTUJ) 
and (llip break the translational and most of the gauge symmetry, leaving only the global gauge symmetry 
given by the gauge transformation 

T 7 : v = (£, a) i-> 7> = (e l7 £, «)■ 



1.7 Stability under gauge periodic perturbations 

We now wish to study the stability of these Abrikosov lattice solutions under a class of perturbations that 
preserve the double-periodicity of the solution. More precisely, we focus on solutions (^,A, $) of ([T]) with 
the initial conditions of the form uq = (^o, A)) = u l + v o> v o € VI ■ 

We now consider the tubular neighbourhood Us of the manifold {T 7 it£ : 7 e 1} of Abrikosov lattices, 
given by 

Us = {T 7 « + v) ■ 7 e R, v e VI, \\v\\ H , < 6} (13) 

We will say that ul is orbitally stable under gauge-periodic perturbations if for all e > 0, there exists 
5 > such that for any C 1 solution A, $), if u = A) e C/5 when t = 0, then it e C/ £ for all i > 0, 
and asymptotically stable, if there exist 7(i) G R, such that any C 1 solution A, $), with u = ('5, A) E Us 
when t — 0, satisfies it — T 7 ( t juJ — > 0, as t — ^ 00. We will say that ul is unstable under gauge-periodic 
perturbations if it is not orbitally stable. 

In order to state are main result, we first need to introduce the Abrikosov function /3(r). We fix b = k 2 
and define /3(t) to be 

^>=wk? (14) 

where £ 7^ is the unique solution of the equation —A A a£ — k 2 £ satisfying the quasiperiodic boundary 
conditions (|10[) . The main result in this paper is the following. 

Theorem 2. For all b sufficiently close to n 2 , the Abrikosov lattice ul is asymptotically stable under gauge- 
periodic perturbations if k 2 > j(l — WF)) an< ^ * s unstable if k 2 < |(1 — jfiyj)- 

To our knowledge there have been no prior results on the asymptotic stability of Abrikosov lattices. 
Orbital stability can often be deduced if the solution in question is a minimizer of an appropriate energy 
functional (see [9]). Hence the orbital stability of Abrikosov lattices for k 2 > | follows from the variational 
proof of [25] that the single vortex Abrikosov lattices for k 2 > i are global minimizers of the Ginzburg- 
Landau energy functional on the fundamental cell (see also |11) . both authors considered only triangular 
and rectangular lattices). However variational techniques do not give asymptotic stability. Also variational 
techniques do not give even orbital stability for k 2 < | as in this case Abrikosov lattices are not global 
minimizers of the Ginzburg-Landau energy functional on the fundamental cell. 

This paper is organized as follows. In the next section we prove the theorem above, modulo a statement 
about properties of the Hessian of the Ginzburg-Landau energy functional. The latter statement is proven 
in Section [3] In the Appendix we use energy methods to prove the weaker statement of orbital stability. 
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2 Proof of Theorem [2] 



As we consider r fixed, from now on we do not display it in the notation. To simplify the notation we 
set 7=1 and cr = 1. To extend the proof of Theorem [2] to 7 and a satisfying Re 7 > and a > is 
straightforward. 

2.1 Decomposition 

The one-dimensional manifold of solutions A4 — {T^Ub :76ft} lies in Us for all S. The tangent space at 
lib € M. is spanned by the infinitesimal global gauge transformation, Tb = 9 7 T 7 M(,| 7= 0j given by 



We now prove the following decomposition for u close to the manifold M.. 

Proposition 3. There exist So > and a map 77 : Us —> K suc/i i/iai T v ^u — Ub _L r& respect to the 

1? inner product). 

Proof. Let X be the affine space X = Ub + Vb- The desired 7 = r\(u) solves the equation / (7, u) = 0, where 
the map / : R x X — > K is defined as 

/(7j u ) = ( T b, T 7 u - u b ) L 2 . 

We wish to apply the Implicit Function Theorem to f(j,u) = 0. It is clear that / is a C 1 map and that 
/(0, Ub) = 0. A calculation shows that D 7 /(0,M{>) = f n \^b\ 2 7^ 0. The Implicit Function Theorem then 
gives us a neighbourhood V of Ub in X and a neighbourhood W of in K and a map /i : V — ?> W such that 
/(7, u) = for (7, zt) £ W x V if and only if 7 = h(u). We can always assume that V is a ball of radius 8q. 

We can now define the map rj on Us for <5 < do as follows. Given u £ Us, choose 7 such that u = T 7 u' 
with tt' € V. We define ?7(u) = h(u') — 7. 

We show that 77 is well defined. We first show that if I7I is sufficiently small, then h(T 7 u) = h{u) — 7. 
First note for all 7, T 7 (V) C V. One can easily verify that f(h(u) — 7, T 7 u) = f(h(u),u) — 0, and therefore 
by uniqueness of H, it suffices to show that h(u) — 7 £ W, but this can easily be done by taking V to be 
smaller if necessary. 

Suppose now that we also have u = T 7 >(ub + v'). Then T 7 '_ 7 (uh + v') — Ub + v. We therefore have 

h(u b + v) - 7 = h(T 7 '-j(ub + «')) - 7 = /i(ufc + w') - (7' - 7) - 7 = /i(tt& + v') - 7', 

so 77 is well-defined. Finally we compute (T b ,T,^ u )U-u b ) L 2 = (T b , T fe ( u /)T_ 7 u— Ub)^ = (Pb, T h {u') u ' ~ u b) = 
f(h(u'),u') = and the proof is complete. □ 

2.2 Hessian 

The chief tool in the proof of the stability result is the analysis of the associated Hessian, i.e., the second 
derivative of the energy functional 8. We first note that £ is a well-defined functional on U§, and that £' is 
in fact its i 2 -gradient in the sense that for all v £ Vb, 



where D is the Gateaux derivative, D£(u)v = d s £(u + sv)\ s —q. We define the Hessian at Ub to be the 
operator Lb := D£'(ub), where D is the Gateaux derivative on maps, on the domain D(Lb) := Vb- Explicitly, 



r 6 = (i* 6 ,o). 



(15) 



D£(u)v 



(£'(u),v) L 2. 




LbV 



( 



Aa 6 £ + K 2 (2|* b | 2 - 1)£ + n 2 ^ b 2 l + 2i{V A i*b) ■ a 
curl* curia + |* b | 2 a - Im(fV Ab * 6 + * 6 V^^) 



) 



(16) 
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This is a real-linear operator on the space L 2 (£l, C) x L 2 (ft,R 2 ) with the domain Vb- We define for it the 
notion of spectrum and of discrete spectrum in the usual way. We introduce the new parameter 



k 2 [(2k 2 - 1)/3(t) + 1] 



(17) 



The term (2k 2 — l)/3(r) + 1 in denominator is necessary in order to have a positive expression under the 
square root and to regulate the size of the perturbation domain. The main result concerning this operator 
which we use in the proof of Theorem [5] is the following theorem, whose proof we postpone until Section [31 

Theorem 4. Suppose that k 2 ^ i(l — ^tj) and that b is sufficiently close to k 2 . Then we have the following 
statements: 

1. The operator Lb has a real, discrete spectrum which includes the eigenvalue with the eigenfunction 
Tb, while its lowest eigenvalue, 9, on the subspace {v G Vb \ v _L Tb} is of the form 

9 := 2b [(2k 2 - 1) /3{t) + l] e 2 + 0(e 3 ). (18) 

Consequently, if k 2 < |(1 — ^^j); then Lb has a negative eigenvalue. 

2. If k 2 > i(l — -gj^pj), then there is a uniform constant c > 0, such that for all v G Vb satisfying v _L Tb, 

{v,Lbv) L 2>c9\\v\\ 2 m . (19) 

3. There exists a positive constant c > such that, for all v G Vb, 

\(v,L b v) L 2\<c\\v\\ 2 Hl . (20) 



2.3 Asymptotic stability 

We now assume that we have the case k 2 > |(1 — jjjkpj) and prove the asymptotic stability of the Abrikosov 
lattice Ub- To this end we derive and use differential inequalities for the Lyapunov functional 

A(t) = ~(v,L b v) L *. (21) 

Lemma 5. Let (^(t), A{t), $(£)) be a solution of the Gorkov-Eliashberg-Schmidt equations (TTJ on the time 
interval [0,T] such that u{t) := {^(t),A{t)) satisfies T l(t) u(t) = u T b +v(t), j(t) G R, v G C 1 ([0, T];Vb), 
v(t) _L T 6 and 11^(0)11^1 « 1. Then \\v{t)\\ H i < e-**/ 2 ||w(0)||ffi for S < 9 for all t G [0,T], where 9 is as in 
Theorem [JJ 

Proof. We plug the decomposition u(t) — T_ 7 ( t )(wb + v(t)) into ^ and use the covariance of this equation 
with respect to the transformation T 7 ( t ) to obtain the equation 

8tv = -L b v + N(v), (22) 

where N(v) is the nonlinearity given by 

ly/ \ ( ^Al ■ « + " ■ V A j + H 2 )i + \a\ 2 * b - n 2 (2 Re(%0 - |£| 2 )£ + n 2 \^ b \ 

for v — ( ) . Using equation (f22|) . we obtain 



-9 t A(t) = (v,L 2 v) L 2 - (N(v),L b v) L 2. (24) 
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Using the expression for N(v) and Sobolev embedding theorems we obtain easily the following rough estimate: 
\Mv)\\ L i<\\v\\ m (\\v\\ m + \\v\\ 2 H i), (25) 
which implies that 

\(N(v),L b v) L *\ < \\L b v\\ L 4v\\ H 2(\\v\\ Hl + \\v\\ 2 Hl ). (26) 
Since, by (|20|) . IHI/p < ||L b t;|| L 2 + ||w||l 2 and ||u||#i < A(t), the last inequality implies 

\(N(v),L b v) L 2 \ < (\\L b vf L2 + A(t))(A(t)* + A(t)). (27) 

Next, we think of L b as the restriction to the orthogonal complement of T b and define L b , for < a < 1, by the 

formula = C J Q (A — L ^ )uj a doj, where C _1 := J Q (^ — Y^j)w a cZu;. (Another way to proceed is to use 

i i 

the complex-linear extension of L b constructed in Subsection l3.3h Then writing (w, L\v)u2 — (L^ v, L b L^ v) ^ 



b u / L - ",^ b ^ b 

and using (fl9| . we find (w,L^u)^2 > 9(v,L b v)L2 = 29A(t). Using this, we obtain 



-d t A(t) > 6A(t) + [i - c(A(i)3 + A(i))] \\L b v\\l, - A(t)i - A 2 (i). (28) 

If we now assume that A(t) <§C 1, then this gives 

-d t (e St A(t)) >(9- S)e St A(t). (29) 
Integrating the last inequality from to t, one finds 

A(0) > e 5t A(t) + (9-5) f e Ss A(s)ds (30) 



and in particular A(t) < e~ 5t A(0) for 5 < 9. Taking A(0)) < 1, we see that our assumption A(t)) 1 is jus- 
tified, which completes the argument. Finally appealing again to (|2"0"|) shows that < e~ 5t ' 2 ||v(0)||#i 
for S < 9. □ 

To finish the proof of asymptotic stability, let So be given by Proposition [3] and let u(0) € By 
standard parabolic existence theory (see e.g. [2TJH2]) for the equation (|22l) . written for the real and imaginary 
parts of v, has a unique solution u(t) € Ug a for t < T, for some T > 0. Let T* be the supremum of such T. 
If < oo, then uiT*) g 9C/,5 . Then by Proposition [31 there is j(t) € R so that T 7 ( t )u(i) = u& + v(i), with 
u(t) € CHtOjT];^) and v(t) _L r 6 . By Lemma H < e- w / 2 ||u(0)|| H i < \5 for all i € [0,TJ, which 

contradicts u(T«) € dE/a . Hence T* = oo and e f/« /2 for all t and ||w(*)||hi < e~ 5t / 2 \\v(0)\\ H i -> 0, as 
< — » oo. Since w(i) = T_ 7 ^u b + T ^{^(t), we obtain (after replacing —j(t) by j(t)) 

\\u(t) - T l(t) u b \\ H i < e~ st / 2 \\u(0) - T l(t) u b \\ H i -> 0, (31) 
as i — > oo, and this completes the proof of asymptotic stability. □ 

2.4 Proof of instability 

We now assume that we have the case k 2 < ^(1 — yjTfj) an d prove the instability of the Abrikosov lattice u b . 
Let r\ £ V b be the negative eigenvector of Lft, corresponding to the eigenvalue —A = 9 < 0, given in Theorem 
IH so that L b rj = — Xr]. We normalize r\ so that 1 1 77 1 1 ^2 = 1. 

For 5 > we now define u{t) to be the solution with the initial datum us.o = u b + Srj. We write this 
solution as u(t) — u b + v$(t). Then vs(t) satisfies the equation (I22j) with the initial condition vs,o = Sr/. 
Using the Duhamel principle and the fact that e~ Lbt rj = e~ xt r/ we rewrite the latter equation in the form 

v s (t)=Se xt r)+ [ e- Lb ^-^N(vs(s)) ds, (32) 



n 
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where N(v) is the nonlinearity given in ([25)1 . It satisfies the following estimate: 

\\N(v)\\v < \\v\\ 2 Hl + \\vf Hl . (33) 

Next, for an appropriate large constant c > A, we have by the standard elliptic theory, similarly to (|20l) . 
IIHIj^i — ( w i(Lb + c ) w )l 2 - The self-adjoint complex-linear extension of the operator Lb obtained in the 
next section and the invariance of the image of Vb under this extension imply the spectral decomposition for 
Lb, which gives that (w, (Lb + c)e~ 2Lbt w) < (—A + c)e~ 2At ||w|| 2 2 . The last two estimates imply the bound 
\\e" Lbt w\\H 1 < e A *||w|j i 2. Using the latter bound and ([33]) and writing in the rest of the proof \\v\\ for \\v Wh 1 , 
we obtain 

\\v s (t) - 6e xt r)\\ < [ \\e- L "^N(v s (s))\\ds (34) 
Jo 

< [ e x ^\\N(v s (s))\\ L ,ds (35) 
Jo 

< ( e^- s H\\v s ( S )\\ 2 + \\v 5 (s)\\ 3 )ds. (36) 



jo 

Now, let M be a constant satisfying < M < \\rj\\ = 1 and define 

Ti := sup{s : \\v s (s) - 5e Xs r)\\ < MSe Xs }. (37) 
Clearly, T\ > 0, and so for < t < T\, we have by triangle inequality 

IK0)|| < \\vs{s) - Se^^W + Se^M < (M + l)Se Xs . (38) 
Therefore, by (jMJ) and (J3HJ) , we have for < t < T 1} 

\\v s (t) - Se xt V \\ < [ e A (*- s )(<5 2 e 2As + S 3 e 3Xs )ds 
Jo 

= S 2 e xt f e Xs (l + Se Xs )ds. 
Jo 

We choose T 2 to satisfy Se XT2 = 1. Then 

\\v 5 (t) - 5e xt r]\\ < CS 2 e 2Xt , for < t < min(Ti,T 2 ). (39) 
Pick C > M, C and define the constant T s > by the relation 

C'5e XTS = M. (40) 

Note that, since C > M, we have T 5 < T 2 . We claim that T s < T x . If not, then T x < T s < T 2 , and by (JUJ) 
and POT) . 

WvsiT^-Se^^W < MSe XTl . 

But this result contradicts the definition of T x in (07]). Hence T 5 < T x , and therefore T s < min(Ti, T 2 ). Now 
we have by the triangle inequality, (IMl) . (|4U)) and the condition 1 1 77 1 1 = 1, 

IMT 5 )|| > WSe^rtW - \\v s (T s ) - <5e Ar %|| 
> Se XTS -6e XT5 M. 

If we set v := (1 — M)^ > 0, then the last equation, together with (14TJ1) . implies ((^(T 5 )!) > which can 
be rewritten as 

\\u(T s ) - u b \\ > v. (41) 

Now we note that for S sufficiently small, the unique minimizer 7* of \\u — T 7 Ub\\ 2 satisfies 7* = O(S) and 
therefore (|4"Tj) implies 

iai\\u(T s )-T^u b \\>u-0(5). (42) 

7 

In other words, V(5 > sufficiently small, there is uq 6 Us such that u(T s ) ^ Ui v , Vt > 0, for a fixed v > 
independent of 5. This implies instability. 
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3 Estimates on Hessian 

In this section we prove Theorem [4] concerning the positivity of the Hessian L b . In what follows we omit the 
subindex L 2 in the inner products and norms. 

3.1 Shifted Hessian 

We first consider a shifted Hessian L b , which induces the same quadratic form. Note that since diva = 0, 
we have curl* curia = — Aa. Again using diva = 0, we can integrate by parts to see that 

f -a-Im(£VA b * 6 + *i,VA;0) = / -2a • Im^V^b). 

Therefore we have (v, L^v) L 2 — (v, L^v) i? where L b is the operator on the same domain Vb given by 

= / -A A „£ + K 2 (2|* b | 2 - IK + k 2 * 6 2 £" + 2i(V Ab * b ) ■ a, \ , . 

bV ~{ -Aa+|* b | 2 a-2Im(£V Ab * b ). ) ^ 



It will be more convenient to work with L, 



b- 



3.2 Rescaling 

We now rescale the problem in order to exploit the analytic properties of the Abrikosov lattice solu- 



tions. Given a pair (^,A) of type (r, b), we set a := = ^y^Hr, and introduce the rescaling 

U a : ($f(x),A(x)) i — ^ (a^(ax),aA(ax)). This has the effect that the rescaled state, (ip,a) := U a ($>,A), 
is of type (r, 1). It is easy to verify that U a is a linear unitary bijection between Vb and V\ (in particular, it 
preserves the L 2 inner-product, i.e., (U^v, U a v') — (v,v f ). Moreover it preserves the orthogonality condition 
in the sense that for v € Vb, v _L r b if and only if U a v _L Ti). 

We note that the rescaled Abrikosov lattice solution (?/>&, a b ) := U<j(^b, Ab) satishes the rescaled Ginzburg- 
Landau equations 

((-A a -x b )^+K 2 m 2 ^=o, (44) 

1 curl* curia — Im(ip'V a ip) — 0, 

where A b = and the quasiperiodic boundary conditions 

4>b{x + v) = ei u ' Jx ip b {x), and a b (x + v) = a b (x) + - Jv, (45) 

where v is either of the basis vectors of C\, as well as diva b = and (a b — ^Jx)d — 0. Thus (tpb, a>b) are of 
type (r, 1). 

We now define the rescaled Hessian to be L r b csc := <r 2 l) a LblJ "o- 1 • With v = ( ^ ^ , it is explicitly given 

by 



n 



r re S c„ _ ( ~A,»£ - A b £ + 2 K 2 |V>f,| 2 £ + «W£ + 2ia • V Q >6 \ 

b \ -Aa + |^ b | 2 a-2Im(£V Q6 ^) )' K ' 

For the rest of this section we write C, ft, V for C\, fii, V\. 

3.3 Complexification 

In order to freely use the spectral theory, it is convenient to pass from the real-linear operator L b csc to a 
complex-linear one. To this end we complexify the space V and extend the operator L r b esc to the new spaces. 
We first identify a : K 2 — >• R 2 with the function a c = a\ — ioti : K 2 — > C. (Whenever it does not cause 
confusion we drop the c superscript from the notation.) We note that a • a' = Re(a c a' ). We also introduce 
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the differential operator d — d Xl ~ id X2 . We note that Ba c = diva — i curl a. where the 8 denotes complex 
conjugate operator. In general, for an operator A, we write A := CAC, where C denotes complex conjugation. 

We now consider the complex Hilbert space L 2 (r2,C 4 ) of vectors (£, <p t a,cv), with the usual L 2 inner 
product 



(v, v') = (i? + #' + aa' + Quj') 



c 



The original space L 2 (ft, C) x L 2 (fl, R 2 ), on which L T b csc is defined, is embedded in L 2 (f2, C 4 ) via the injections 
7r±: {Z,a)^^=(Z,±i,a,±a), (47) 

with inverses, tt^. 1 , given by the obvious projection. V± := Ran7r± are real spaces spanning L 2 (R 2 ,C 4 ): 

(£,<p,a,u) = + + (j),a + uj,a + uj) + ^{£_- + <j),a- Q,-a + u). (48) 

Moreover, the map I : (£, <fr, a, uS) — ¥ (i£, icj>, ia, iuS) acts between V±: I : V± — > V^. This embedding transfers 
the operator U b esc to the range, V± := Ran7r±, of this injection. We denote the resulting operator by L)j ransf . 
Its domain is ttV '. We want to extend it to L 2 (R 2 , C 4 ). To this end it is convenient to rewrite the operator 
£transf - n com pi ex notation. We introduce the notation d a c = d — ia c . Straightforward calculations show 
that 

2ia ■ V ab ip b = -i{d* a c4i b )a c + i(d a c^ b )a c , 
and that 



Using the above relations we rewrite the operator L^ ransf and then define its complex-linear extension, denoted 
by K b , by the resulting matrix 



K h = 



/ -\ b -X b + 2K 2 \^ b \ 2 _ K 2 ip 2 i( d a M \ 

* 2 j>l -/S ab -X b + 2n 2 \i, b \ 2 ~i(d a M i(d*J*) 

i{d a M -A+|^ h | 2 

V -%(Mb) -i{dl b i>b) -A + |V> 6 | 2 / 



(49) 



on the domain which consists of all v = (£,<p,a, uj) € _ff 2 (M 2 ,C 4 ), with £, <j>, and a, u> satisfying the quasi- 
periodic boundary conditions 

X(x + v) = e5 l " Jx ^(x), and a(x + v) = a(x), (50) 

where, as above, v is either of the basis vectors of £, as well as a and w are divergence free and have 
mean-zero, diver = and {<r)c = 0. (Note that similarly to the Riesz - Fischer L 2 — space on a torus (see e.g. 
[21]). we could have used results of [35] to introduce L 2 — space on fl satisfying the quasiperiodic conditions 
([50)) . rather than periodic ones.) 

The operator K b is clearly complex-linear, self-adjoint, has purely discrete spectrum, and, as it is not 
hard to check, satisfies 

V± are invariant under K b and K b \ v+ = L£ ransf , (51) 

K b = K b \ v+ +K b \ v _ 1 a(K b \ v+ ) = a{K b \ v _) and a(K b ) = a(K b \ v+ ) U a(K b \ v _). (52) 
For the second equation, we used (|48p and that K b obviously commutes with I. (|5 1 [) implies 

(v, Lr c v) = (<k+v, Lf^n+v) = (n + v, K b n+v) . (53) 
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3.4 Perturbation theory 

It is shown in [32 that for each r there is eo > 0, such that the solutions (ipb, a b) form a real-analytic branch 
of solutions in the (bifurcation) parameter e € [0, eo] defined in (|17j) and have the following expansions for 
(see [32] where a more general result was derived): 

ip b = eil>° + e 3 ip 1 +0(e 5 ), 
a b = a° + e 2 a 1 + 0(e 4 ), 

where ip° and a 1 satisfy the following relations 

«°=0, <|V>Tk=2, (54) 



l da} = -((\i; \ 2 )c 



Aa 1 = ~<4>\d a oi> ). 



(55) 
(56) 



We now prove the following proposition. 



Proposition 6. Let be the gauge zero mode in the extended space: T£ = (iipb-, —iipb, 0, 0) € V+ . For e > 



sufficiently small, the lowest eigenvalue, 9, of K-\ 
is of the form 

9 := 2 [{2k 2 - 1) /3(r) + l] e 2 + 0(e 3 ). 



A" 



trans f 



on the subspace {v €E D(K+) \ v _L T} 



(57) 



is afeo i/ie lowest eigenvalue of L r b esc on (rjj) and therefore b9 is the lowest eigenvalue of L b on 



Proof. We use the expansions above to expand Kf, in powers of e and the relation d* ^° = 0, to simplify the 
resulting terms to obtain K b = K° + eW 1 + e 2 W 2 + o(e 3 ), where 



( -\o - 1 




V 



w = 














-A^-l 
















\ 










-A 










-A 





(58) 



KdaOlp ) \ 

-*(cW°) o 



w = 



K 2 (^ ) 2 






K 2 (<) 2 

B° 



















IVT J 



J 



where 



-A 1 + 2k 2 \^\ 2 - ia 1 d* a o + i^dao. 



(59) 



(60) 



(61) 



The unperturbed operator K° reduces to the operators studied previously, e.g. in [32], where these 
operators are denoted by L n with n = 1 and M and where the spectra of the latter operators are described 
in details. In particular it is shown there that, in general, K° has the eigenvectors (V>°, V' , 0, 0), -0°, 0, 0), 
(0,0,1,0) and (0,0,0,1). The latter two are ruled out by the condition that (a) = 0. We summarize the 
properties of K° in the following lemma. 
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Lemma 7. K° is a nonnegative self-adjoint operator with discrete spectrum. It has a zero eigenvalue of 
multiplicity 2 and the kernel is spanned by the elements 

« 0p = W-°,/^ ,0,0), p = ±. (62) 

Note that vq± € V±. Hence the operator K+ := K°\y + which is the zero-order approximation of the 
operator K+ := Kb\v+ that we are interested in, has the simple lowest eigenvalue with the cigenfuction 
v a +. Also, the zero mode vq- is related to the complexified gauge zero mode := n-Tb — (iipb, —i^Pb, 0, 0). 
Indeed, expanding the latter vector in e, we obtain 

T c b = iev Q . + 0(e 3 ). 

By standard perturbation theory (see e.g. [551 [TSl E3), the spectrum of K b consists of eigenvalues 
which cluster in e— neighbourhoods of the eigenvalues of K° and each cluster has the same total multiplicity 
as the eigenvalue of K° it originates from. Of course, the same is true for its restriction K+ := Kb\v+- 
Namely, the spectrum of K+ consists of eigenvalues which cluster in e— neighbourhoods of the eigenvalues 
of K+ := K°\v + and each cluster has the same total multiplicity as the corresponding eigenvalue of K+. 
Thus K+ has the simple eigenvector Vb, which is a perturbation of the simple eigenvector vq+ of K+ with 
the smallest eigenvalue 0. It suffices to determine the corresponding eigenvalue of K+, which we denote 6. 
Since Kb is self-adjoint and rg is its eigenfunction with eigenvalue 0, we have that Vb -L T. (In the leading 
order this becomes i>o+ _L Uq_ which can be verified directly.) Hence is the smallest eigenvalue of K+ on 
the subspace orthogonal to T%. 

To find 9 we use the Feshbach-Schur map argument (see e.g. [16[ 17]) with the projection P given by the 
orthogonal projection onto vo+ G null if?.. This argument implies that 

A <= o-{K+) if and only if A e cr(J>(A)), (63) 

where, with P = 1 — P. 

7>(A) := [PK+P - PK + P(PK+P - X^PK+P]^ , (64) 

provided the operator PK+P — A is invertible on RanP and PK+P is bounded. (The latter conditions 
suffice for the right hand side of (|64[) to be well defined.) Due to the relation PKbP — P(K+ — K+)P and 
the straightforward estimate 

\\K+-K°+\\<e, (65) 

we see that the operator PK+P is bounded. To show the invertibility of the operator PK+P — A on RanP, 
we note that it is the restriction of the operator QK+Q — X, where Q is the orthogonal projection onto null K° 
and Q := 1 — Q, to the subspace V+ We know that a(QK Q) C [vq, oo) for some Vq > and therefore, by 
standard perturbation theory we have that 

a{QK b Q\ Rlin Q) c[c,oo), (66) 

with c = Uq + 0(e). Hence the self-adjoint operator QKbQ — X is invertible on RanQ, provided A < c, and 
therefore its restriction PK+P — X (to real-linear subspace V+) is invertible on RanP and \\(PK+P — A) _1 | < 
c _1 (again provided A < c). Hence (fB"4")) is well defined for A < c. 

We now use K ± = K + +eW 1 +e 2 W 2 +o(e 3 ), the relation if°P = PK% = and the facts ||PJf+P|| = O(e) 
(by fBjJ)) and || (PK+P — A) _1 j| < 1, provided A < c (by Since we are studying the eigenvalue in 

0(e)— neighbourhood of 0, we have that A = O(e). Using this, we obtain 

J"(A) =eJ- 1 +e 2 J- 2 + 0(e 3 ), (67) 

where 

T x := (i; 0+ ,H /1 «o+)/(|i'o+| 2 >£, T% ■= («o+, [W 2 — W 1 P(PK+P)^ 1 PW 1 ]v +) / (\v +\ 2 ) C - (68) 
The operator, W\ is explicitly given by ([551) . This expression implies that (vq+, WiVq+) = and therefore 
J 7 ! '■= (vq+ 1 W 1 v +) = 0. (69) 



We now turn to the e 2 order operator, Ti. 
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Lemma 8. 

F 2 = 2[{2k 2 -\)(3{r) + l]. (70) 



Proof. We begin with (vq+, W 2 v +). We first note that W 2 and vq p are explicitly given by (JHOJ) and (f62]) . 
Using the fact that d* ip° = 0, we calculate that 

(v 0p , W 2 v 0p ) = 2 (-A 1 (|^°| 2 ) £ + 2k 2 (\i/j°\ 4 ) c + Re^ia^V + P k 2 (\^ \ 4 )c) , (71) 

where Ax is defined by the expansion A& = 1 + e 2 A x + 0(e ) and, due to the definition of Af, and e, is equal to 

A J = 



^ + (* 2 -^)/3(t) 



(IV>T>£. (72) 



We note that —A = 8*d. Using the identities (|56t and (|55|) we obtain 

2Re( l a 1 ^°cWV = 4(a 1 Aa 1 ) £ = -4(|cV| 2 ) £ = -(|V>T>£ + (l^°| 2 ) £ - (73) 
The equations (fTTj) and (|73]l and the relation ffT2]) give 

(v 0p ,W 2 v Qp ) = 2(l + p) K 2 (\^ \ 4 ) c . (74) 
To compute the second term in T% we note that PW 1 P = W 1 P ', and use (f58|) . we calculate 



(uop.W^Pi^PrW^op) = -2Re(^ (a a oV )A- 1 (^ (9 a o^))) £ . (75) 
Again using the identities (f56|) and (j55|) and the fact that —A = 8*8, we obtain 



-2Re(V; o (a a o^)A- 1 ^ (9 a oV )))£ = -8Re(a 1 Aa 1 ) £ = S^a 1 ! 2 ) 



= 2(|^| 4 ) £ -2(|^| 2 ) 2 . 

The second equation in (|54|) and the fact that (|wo P | 2 )z; = 2(|V' | 2 )z: now gives (f70|) . □ 



(76) 



The equations (|63)l . (j67|) . (|69| and (j70|) imply the first part of Preposition [6] Finally, by the definition, 
9 is also the lowest eigenvalue of L r b csc on T] 1 and therefore, by the formula relating L r b esc to Lb, we see that 

b 



Ob is the smallest eigenvalue of Lb on . □ 



Now we are ready for 

Proof of Theorem^ We restore the subindex L 2 in the inner products and norms. Preposition [5] implies the 
first statement of the theorem and the estimate 

(v,K b v) L 2 > e\\v\\ 2 L2 . (77) 

for all v £ D(Kb), such that v _L Tg. The latter estimate gives, for all v £ D(Lb), such that v _L Tb, 

{v,L b v) L .> l -0\\v\\ 2 L2 . (78) 

We upgrade now the lower bound on (v, LbV)i,i to that on (v, Lbv) #i . (We could have done with the operator 
Kb as well.) We begin with 

Lemma 9. For all v £ Vb, if v _L Tb, then 

\\\v\\ 2 m-C\\v\\%<{v,L b v) L *<\\v\\ 2 m , (79) 
for some positive constant C . 
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Proof. We write v — (£, a) and also A J = A® + P. For convenience we simplify the notation in the following 
calculations. Integrating by parts we obtain 

(v,L b v) L * = -L Re ^ -£A A o£ + 2iP • £V A o£ + |P| 2 |£| 2 - * 2 |£| 2 + (2k 2 + i)|* b | 2 |£| 2 



1. 

2a ■ Im(£V A o W) - 2(a • P) Re(£* 6 ) 



+ (« 2 - + 2l « • (£V4?*b) + 2(a ■ P)* b £ -a-Aa + |* b | 2 |a| 2 



= P" 1 IVa ^ + |Va ' 2 + |P|2|e|2 " K2|?|2 + ^ + ^> )|VN2|e|2 



la 

Using this expression and the estimate obtained with the help of the Schwarz inequality 



+ (« 2 - 5)*?? + l*6| 2 M 2 - ^ • hniLY x i i 4o -Tn.liV, (S(i, 



1 P-Im(^oO 
for any r > 0, we obtain 



< ||P||oc||?|UH|V A oe|U= < JllPHoo (lIKI^ +r||V A o£||£ 2 ) , (Nli 



-?lieil! 2 -rC||V j4 oC||| 2 -C||a||| 2 -C||e|| 



£2- 



Now we choose r as r = ^ so that we arrive at the lower bound in (|79l) . To obtain the upper bound we use 
(|50)l and (EJ) again and the fact that ||*b||oo, ||P||oo < oo. □ 

Let now S £ [0,1] be arbitrary. Then using (|77|) and (|79|) 

(v,L b v) L 2 = (1 - 5)(v,L b v) L 2 +S(v,L b v) L 2 

>(l-5)6\\v\\l 2 +5(±\\v\\ 2 Hl -C\\v\\l 2 ) 

= ((l-5)e-dC)\\vf L2 + 5 -\\vf Hl . 

(IT9l) now follows by choosing <5 = i + g +c - Next, the estimate (j20|) follows from the upper bound in ([79| . □ 

A Proof of orbital stability 

We now assume that we have the case k 2 > |(1 — ^tj) and prove the orbital stability of the Abrikosov 
lattice u\. This is a weaker statement than the asymptotic stability which we have already proven but it 
requires rougher analysis. We follow [14] . 

Theorem 10. For all b sufficiently close to k 2 , the Abrikosov lattice u r h is orbitally stable under gauge- 
periodic perturbations if k 2 > ^(1 — WF))- 

Proof. As in the main text, we consider r fixed and do not display it in the notation. We will require a series 
of lemmas. 

Lemma 11. There exists positive constants c and C, such that for all v £ V b , if v _L r b , then for any 9' < 6, 
0>\\m - 4vfm - c\Hm < ^(u b + v)- £{u b ) < C{\\vf m + \\v\\ 3 H1 + \\vf H1 ). (82) 
Proof. Using the Taylor expansion of £, together with the fact that £' (u b ) = and v _L r b , we have 

£(u b + v) - £{u b ) = ~(v, L b v) + R(v), 
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where the remainder R(v) is given by, setting v = (£, a), 

R(v) = / M 2 (Re(* b £) + h^ 2 ) ~ a ■ Im(fr Ab O + ^(M^bO + ^| 2 )- 

Using the Cauchy-Schwarz and Sobolev inequalities it is straightforward to show < c(|| t'll^i + 

Using (fT!)|) . this gives 

£{u b + v) -£{u b ) > 9\\v\\ 2 Hl -c\\vf H i -c||u||^i. 

On the other hand, by the definition (|16p of L b and the boundedness of u b together with its derivatives, 
we have (v,L b v) < This estimate together with the above estimate of gives the upper bound 

in and this completes the proof. □ 

Lemma 12. Suppose that (ty, A, $) is a solution of the Gorkov-Eliashberg- Schmidt equations ([T]) on the time 
interval [0,T] satisfying (I) - (III) and u(t) = (^(t),A(t)) £ C^^T]; Us) for any S > and T > 0. Then 
the energy function £(u) is a nonincreasing function in time. 

Proof. Note that we have dtu £ V b . Using the gradient-flow form of the equations (|4]), we see that 

d t £ (u) = (£'(u),d t u) L 2 

= -(£'(u),(iW,V$)) L *-\\£'(u)f L > 

< I Im($^A A *)-Re(iK 2 (l-|*| 2 )|*| 2 $)-curl*curlA-V$ + Im(§Vi*)-V$ 

jQ b 

= I -Im(WA*)-V$-Im($|VA*| 2 )+Im(*VA*)-V$ 
= 0, 

where we used the fact that $ is real- valued and div curl* =0. □ 

Lemma 13. Let A(t), $(t)) be a solution of the Gorkov-Eliashberg- Schmidt equations ([lj on the time 

interval [0,T] and denote u(t) — (^(i), A(t)). Given e > sufficiently small, there exists 5 > such that if 
T l(t) u(t) = u b + v(t), 7 (f) e R, v £ C-dO^^Vh), v(t) _L r 6 , and \\v(0)\\ H i < 5, then \\v(t)\\ H i < e for all 

t e [o,T]. 

Proof. We set N(t) = ||i>(i)||fl-i. Using the inequalities (|82|) and Lemma H2l we have 

C x N{tf - C 2 N{tf - C 3 N(t) A < £(u b + v(t)) - £{u T b ) = £ (T 7(t) u(t)) - £(u b ) 

= £{u{t)) - £{u b ) < £(u(0)) - £{u b ) 
= £ (T 7(o) u(0)) - £{u b ) = £{u b + v(0)) - £{u b ) 
< C 4 N(0) 2 + C 5 N(0) 3 + C 6 N(0) 4 . 

Now there exists S > such that if the left-hand-side CiN(t) 2 - C 2 N(t) 3 - C 3 N(t) 4 < S , then either 
< N(t) < e or N(t) > e' for some e' > e. We can choose 6 sufficiently small so that the right-hand-side 
C 4 jV(0) 2 + C 5 iV(0) 3 + C 6 7V(0) 4 < S . The result then follows from the continuity of N(t). □ 

We can now prove the following proposition, which implies Theorem 1101 

Proposition 14. For any e > 0, there exists S > such that if (ty, A, $) is a C 1 solution of the Gorkov- 
Eliashberg- Schmidt equations (JTJ and u(t) = (\t (t) , Ait)) satisfies u(0) £ Us, then u(t) £ U e for all t>0. 

Proof. Let <5o be given by Proposition [3] If u(Q) £ Ui Sq , there is T > such that u(t) £ Us for t <T. Then 
by Proposition H there is y(t) £ 1 so that T l(t) u{t) = u b + v(t), with v(t) £ C 1 ([0, T];V b ) and v(t) _L T b . 
By Lemma [TBI we can then find Si such that ||u(i)|| < ^min(e,<5o) for all t £ [0,T]. 

Now let S = min(5o,^i). It follows that if u(0) £ Ui s , then u(t) cannot leave U e D Us , and that proves 
the proposition. □ 

This finishes the proof of Theorem [TUJ □ 
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